% This m-file is to export the isotherms
close all
clear all
clc
IsothermsFile={
    'sph01_scan106_iso.txt';   % #1 KCl w/o Fe
    'sph02_scan21_iso.txt';  % #5   KCl/FeCl3 pH 2
    'sph01_scan171_iso.txt';   % #2 KCl/FeCl3 pH 3
    'sph01_scan315_iso.txt';  % #4  KCl/FeCl3 pH 4
    'sph02_scan86_iso.txt'  % #6    KCl/FeCl3 pH 5
    'sph01_scan245_iso.txt';  % #3  KCl/FeCl3 pH 7.1
    }
Colors={'k','b','y',[0 0.4 0],'m','r'}
figure
hold on
factor=1/4.31;
for k=1:length(IsothermsFile)
    data=dlmread(IsothermsFile{k});
    area=factor*data(:,1);
    press=data(:,2);
    plot(area, press,'-','LineWidth',3,'Color',Colors{k});
end
hold off
set(gca,'LineWidth',3,'Box','on');
set(gca,'ylim',[-1 34],'xlim',[39,64],'xtick',[40:5:65])
set(gca,'FontSize',30);
axis square